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Abstract 

We consider the higher order gravity with dilaton and with the leading string the- 
ory corrections taken into account. The domain wall type solutions are investigated 
for arbitrary number of space-time dimensions. The explicit formulae for the fixed 
points and asymptotic behavior of generic solutions are given. We analyze and classify 
solutions with finite effective gravitational constant. There is a class of such solutions 
which have no singularities. We discuss in detail the relation between fine tuning and 
self tuning and clarify in which sense our solutions are fine-tuning free. The stability 
of such solutions is also discussed. 
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1 Introduction 



Theory of closed strings incorporates gravity in a natural way since the spin-2 field is always 
present in its massless spectrum. From this point of view theory of gravity is an effective 
theory of the graviton when the energies are smaller than the mass of the string first excited 
level. String theory contains many other massless fields and in an effective action one 
expects the presence of all powers of the fields - indeed the calculations [1] of the string 
amplitudes involving the graviton, dilaton and the antisymmetric tensor confirm that the 
string effective action contains higher order terms. The expansion parameter in string theory 
a' is dimensionful and connected to the fundamental string length l s by a' = l~ 2 . Since the 
effective action was calculated in [1] (to order a' 2 ) by comparison with the string amplitudes 
on-shell, we cannot distinguish in this way actions differing by lower order equations of 
motion. However these actions from the field theory point of view are different and the 
equations of motion have different solutions. Since we don't know string amplitudes off-shell 
we have to employ other methods to try to guess the true effective action - one possible 
idea is to look for an action with some additional, string motivated symmetry. In [2] it 
was proven that one of the actions for gravity coupled to the dilaton and the antisymmetric 
tensor at order a' 2 exhibits an 0(d,d) symmetry characteristic for string effective actions 
[3]. Therefore in the present paper we use the form of the action obtained in [2] and analyze 
the presence of the domain wall type solutions of the equations of motion. 

It was shown in [4] that in the presence of arbitrary number of Euler densities in the 
Lagrangian (without the dilaton), there always exists a domain wall solution of the Randall- 
Sundrum type [5] . In the present paper we analyze the gravity theory coupled to the dilaton 
with interaction terms up to the fourth order in the derivatives. Properties of the generic 
domain wall type solutions are investigated. We are interested mainly in those solutions 
which allow for localized gravity with finite effective gravitational constant after integrating 
over the coordinate perpendicular to the brane. We find that both with and without the 
cosmological constant the presence of the dilaton allows for solutions of the Randall-Sundrum 
type [5] as well as solutions with singularities of the type described in [6] . We find also some 
new cases specific to higher order gravity. It turns out that some of the solutions without 
singularities have finite effective gravitational constant without fine tuning of the parameters 
of the Lagrangian. The problem of the absence of fine tuning and potential self tuning (for 
the discussion in the case of the lowest order gravity see [7]) of these solutions is addressed 
in great detail. We explain the relation between our fine-tuning free solutions and some 
no-go theorems [8, 9] derived for theories similar to the one consider in the present paper. 
Recently one similar solution for a 5-dimensional theory has been found in ref. [10]. However 
the higher order terms used in [10] are as considered in ref. [11] and differ from those adopted 
in the preset work (we explain the reason for this difference in the next section). Moreover, 
we present more solutions of this type and our analysis can be applied to theories in arbitrary 
dimensions. We discuss also, contrary to [10] , the stability of such fine tuning free solutions. 
Some exact solutions in higher order gravity have been recently discussed in [12]. 
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The paper is organized as follows: In section 2 we present the Lagrangian and derive 
the equations of motion in general and in the application to the brane metric. In section 
3 we explicitly find fixed points and asymptotic behavior of solutions. The number of the 
fixed points and the behavior of asymptotics depends very much on the value of the bulk 
cosmological constant. In section 4 we analyze and classify solutions with finite effective 
gravitational constant - the finiteness is due either to the compactness of the transverse 
dimension or to the decreasing warp factor like in the Randall-Sundrum scenario. Section 5 
is devoted to the detailed discussion of the equations of motion at the brane. In section 6 we 
discuss the problem of fine tuning of the solutions with finite effective gravitational constant 
and clarify in which sense our solutions avoid fine tuning. The relation of our results to some 
no-go theorems [8, 9] is also explained. In section 7 we present the conclusions. 

2 Higher order action and equations of motion 

In theories with gravity coupled to the dilaton one has to decide upon the so called frame (i.e. 
choice of metric) since one can redefine the metric by a (non-vanishing) factor depending 
on the dilaton. There are usually two frames used: the string frame (with the factor in 
front of the curvature scalar in the string tree level effective action) and the Einstein frame 
(without this factor). The metrics solving the equations of motion can have quite different 
properties in the two frames (for a discussion in the cosmological context see [13]). In the 
present Universe and with the presently available energies the difference between the frames 
is unimportant since the dilaton seems to be (by a still unknown mechanism) constant both 
in time and in space, but both in the early Universe or, as discussed in this paper, close to 
a brane, the difference can be quite substantial. It is generally argued that although the 
physical conclusions should not depend on the redefinition of fields, it is much easier to draw 
these conclusions in the string frame than in the Einstein frame and the "translation" from 
one to the other is not always evident. Therefore we will use the string frame in this paper. 
It should be pointed out that the higher order terms considered in refs. [10, 12] are different 
from those used in the present paper even after changing the frame. 

The string motivated action we consider reads [2] (note a difference by a factor of 2 in 
normalization of the dilaton 0): 





where a = a'\ and 




(2) 





(3) 
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Lb = ^y/Hviffliv) ■ ( 4 ) 

Rgb is the Gauss-Bonnet term 

Rgb = Rt ivp ° R^upa — ^R^Rfiu + R 2 ■ (5) 

The space-time has arbitrary dimension D. The metric in eq. (4) is the induced metric on 
the brane 

g a p = g^6£6% . (6) 

The indices have the following ranges: v — 1, . . . , D; a,(3= 1, . . . , d; where d = D — 1. 
The coordinate perpendicular to the brane we denote by y = x D . We fixed the position of 
the brane at y — 0. 

Let us now calculate the functional derivative of this action (separately for each term in 
(1)) with respect to variations of the metric and to variations of the dilaton. We define 



Tff = 2« 2 [d D xL I: (7) 

W (i) = 2k 2 ^L^- [d D xLr, (8) 
5(f) J 

for I = 0, 1, B. A rather tedious calculation gives for T^ u 

T$ = Ag^ + R„ v - -g^R + ^g^d^d^ + - (□</>) ^ , (9) 



T^J — AR kX R k ^\ u — 2R^ crK \R u aKX + AR^ K R K U — 2RR llv + -g^R-GB 
+4(V K d x <f))R tMKuX + W^d^R^ + 8(V K ^0)i?% - 4(n0)i^ 
-4(V K <9 A 0)/T V - 2{d K 4>d«4>)Rg^ - 2{V»d v <j>)R + 2(D0)% i „ 
-2<9 K <R0D0^ + 2(D0) V - 2(V K <9 A 0)(V K <9 A 0)^ + 2d K <j)d K <j)V fl d u <t) 



+4(V,A0)(V„<9 K 0) - 4(V / A0)D0 + -(d«<R0) V , (10) 

?J? = -yMSWtfWi/) • ( n ) 

Calculation of the functional derivatives with respect to is straightforward 

W (0) = 2A + d K <pd K (j) -20(f)- R, (12) 

= R 2 GB + ^ (r kX - \Rg KX ) (d K <pd x <p - 2V K 9 A 0) + (d^dy) 2 

+4(D0) 2 - 4(V K <9»(V K <9 A 0) - 4R KX d K( f)d x( f) , (13) 

W (B) = e <V(0)<%) . (14) 
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The equations of motion read 



T (o) + aT (i) + T {B) = o 
W^+aWW + WW =0. 



(15) 
(16) 



Let us note that in the bulk, i.e outside of branes and singularities, we have the equality (up 
to boundary terms) that can be easily checked 



S = 



1 



2k 2 



d D x^ge~ 4 ' \W {0) + aW w 



(17) 



It means that on-shell i.e. on the equations of motion the effective cosmological constant 
has no contribution from the bulk. On physical grounds, however, we are interested in the 
solutions for which the total (i.e. including branes) effective cosmological constant vanishes 
on-shell. Anticipating the results of section 5, for fixed point solutions it happens for one 
form of the potential on the brane, namely V(<f>) = Xe~ 4 i.e. the potential suggested by 
string theory (for NS fields for which the bulk action is considered). For special solutions 
which are not exactly at the fixed points it happens for larger class of brane potentials. Then 
the boundary terms from the bulk cancel the brane contribution and we have 



2k? 



J d D x^e~* [W<M + aWW + W^] 



(18) 



where the integral is over the full space including brane positions. Since it is proportional to 
the equations of motion, it vanishes on-shell so the effective cosmological constant vanishes 
as well. 

We now start with solving the equations in the bulk and later we will include Tj£) and 
W( B \ The line element of the domain wall background is equal to 

ds 2 = e 2A{v) r ]a/3 dx a dx (} + dy 2 . (19) 

The three equations of motion - a (3 and yy components of (15), and (16) - read respectively: 





A + (d— 1)A" - 0" + l -d{d - 1)(A') 2 -(d- l)A'<f>' + i(0') 2 



+a 



2{d - \){d - 2){d - 3) A" (A) 2 - 4(d - l)(d - 2)A'A^' + 2{d - l)A"(<p') 2 



-2(d- l)(d-2)<f>"{A) 2 + 4(d 



'A 



L"/ 



+-d(d -l){d- 2){d - 3)(A') 4 - 2{d - l) 2 {d - 2){A'f4>' 
2 



+{U - A){d - 1){A 



l\2f il\2 



2(d-l)A'( 



k'\3 



(20) 



= A + ±d(d-l)(A) 2 -dA<f ) > + l( ( j ) >) 



/\2 
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+a 



-d(d - l)(d - 2)(d - 3)(A') 4 - 2d(d - l)(d - 2) (A') 



+3d(d - l)(A'Y((j)'y - 2dA 



= A + dA" - 0" + ^d(d + l)(A') 2 - dA'<j>' + i(0') 2 



1 

2 



A4 



(21) 



2d(d - l)(d - 2)A"(A ) - 4d(d - + 2dA"(0 

1 



i/\4 



-2d(d - 1)0 (A ) + 4d<t>"<t>'A' - 20"(0') + -d(d + l)(d - l)(d - 2)(A') 
-2d 2 (d - 1)(A') 3 0' + d(3d - 1)(A') W - 2dA'(0') 3 + ^(0') 4 
These equations have two obvious symmetries 

V ^V + Vo, A-> A + A , 0^0 + 0o, 

and 



y 



-y, A'^-A', 



(22) 



(23) 



(24) 



which can be used to obtain a whole class of solutions starting from each particular one. 

Only two of the equations (20)-(22) are independent in the bulk. One of the combinations 
of those equations is particularly simple. Multiplying the difference of (20) and (22) by 
exp(0 — dA) we get an expression which is a total derivative with respect to y. Integrating 
it we obtain 



e dA-4, | A / [ x + 2a ^ d _ _ 2) A' 2 - 2(d - l)A'(f)' + 0' 2 )] } = const . 



(25) 



This equation will prove to be very useful when analyzing the fixed points of the equations 
of motion. 



3 Fixed points and asymptotic behavior of solutions 

Let us find the fixed points of the equations of motion (20)-(22) (i.e. values of A' and 0' for 
which A" = 0" = 0). We denote 

A' = ^, 0' = ^, (26) 

where ai, b\ are constants. Substituting constant derivatives A' and 0' into eq. (25) we see 
that it can be fulfilled only in two ways: either 0' = dA' (only then the exponential factor of 
(25) is constant) or the expression in the curly bracket in (25) must vanish. It is not difficult 
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to check that only the second possibility is compatible with eqs. (20) and (21). Thus, every 
fixed point of eqs. (20)-(22) must satisfy 



l + 2(d- l)(d - 2)a\ - A(d - l)ai&i + 2b\ 



0. 



(27) 



There are two types of fixed points corresponding to two ways in which the above equation 
can be fulfilled. The first type has vanishing a±. The constant b\ can be then calculated 
from eq. (21). The result reads: 



di — , bi — ±<7 , 



where 



Vl-8aA- 1 



(28) 



(29) 



The number of such fixed points depends on the sign of the bulk cosmological constant: there 
are two for A < 0, one for A = and none for A > 0. The solution for the vanishing bulk 
cosmological constant has a± — b± — and is just the 5-dimensional Minkowski space-time. 

For the fixed points with nonzero a\ the expression in the square bracket in eq. (27) must 
vanish. This gives a relation between a\ and b\ which together with eq. (21) gives the second 
type of the fixed point: 



ai = e a 7, h = e a (d - i)j + e b \ (d - 1)7 2 



1 

2' 



(30) 



where e a , e b = ±1 and 7 can take one of the two possible values, 7+ or 7_, given by 



7± 



1 ± 



4aA(d-2) 



\ 2(d-2) ^ \ 2(d- 1) d-l 
The requirement that the parameters a\ and b\ are real gives two possibilities: 

7 = 7+ for A > A , 



or 



where 



7 = 7- 



A 



for Ai > A > A , 



8a(d-2) ' 
A 1 = A + [4a(d- l)(d- 2)]' 1 



(31) 

(32) 
(33) 

(34) 
(35) 



If solutions for 7 = 7+ (7 = 7_) exist, there are four of them corresponding to four combi- 
nations of signs e a and e& in eq. (30). The solutions for 7 = 7. exist only for a rather small 
range (33) of values of the bulk cosmological constant. 
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Adding all the above types of solutions we can see that the eqs. (20)-(22) have the 
following number of fixed points: 2 for A < A ; 4 for A = A ; 10 for A < A < Ai; 8 for 
A = Ai; 6 for A x < A < 0; 5 for A = 0; 4 for A > 0. 

Let us now analyze the asymptotic behavior of A' and </)'. There are different types of 
behavior related to each of the allowed fixed points. One can find also asymptotic behavior 
for A' and going to infinity. The results follows: 



There is one fixed point with A' = 
y — > ±oo is given by 



b' if A = 0. The asymptotic behavior for 



A' 



1 



y/d\y\ 



0' 



y/d- sgn(y) 



(36) 



For the two fixed points with A' = for A < and y — > ±oo we have 

A' = a 2 exp(-cr|y|/ y/a) , # = -sgn(y)-^= + da 2 ex.p(-a\y\/y/a) , (37) 

y/a 

where a is given by eq. (29). 

If the bulk cosmological constant satisfies the condition Ao < A there are four fixed 
points A' = ± r y/y/a with 7 = 7+ corresponding to four combinations of signs e a , 
in the expression for b\ given in eq. (30). In addition if A < A < Ai there are four 
fixed points with 7 = 7.. Half of these fixed points are approached for y — > +00 and 
the other half for y — > — 00. Let us analyze solutions approaching the fixed points for 
y — > 00. The solutions can approach them either from above (03 > 0) or from below 
(a 3 < 0) and their asymptotic behavior is given by 
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A' = e a —= + a 3 exp(-cy/ y/a) , </)' 
y/a 



where 



c = e a 7 - e b ]l (d - 1)7 2 - ^ ' 



6 3 /a 3 = d 



1 - 



2(d-l) 1 \e a e b \l(d-l) 1 *--+ 1 



-^= + b 3 exp(-cy/y/^), (38) 
\ a 



(39) 
(40) 



From the condition c > and the fact (cf. (31)) that for d > 2 



{d-l) 1 l--< 1 .. 



(41) 



we conclude that solutions with 7 = 7+ approach (as y — > 00) the fixed point when 
€b = — 1 while those with 7 = 7. when e a — 1. 
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• There are vertical asymptotes for each value of y. Denoting the location of the vertical 
asymptote as yo one can check that \A'\ and \<f>'\ go to infinity for y — > y as 

A' = — ^— , 0' = , (42) 

y - yo y-yo 

where c, are four solutions (two positive and two negative) of the equation 

d(d - 2)c\ + 8dc 3 i - ISdcj + 27 = . (43) 

A 1 is positive for positive (negative) q and y > yo (y < yo). The solutions with negative 
A' are as usual related to those with A' > by the symmetry (24). Let us note that 
Ci < 3/2 for d > 2. 



A 



Y 








y 

Figure 1: Typical solutions of the equations of motion (20)-(22) for d < 9 and the bulk 
cosmological constant in the range A x < A < 0. Only solutions with positive A' are shown. 
The cases with negative A' can be obtained by means of symmetry (24). From the symmetry 
(23) it follows that any of the solutions gives another solution when shifted by arbitrary 
constant along the y coordinate. The circles denote singularities of type B (see discussion in 
section 4) for which the metric is non zero but the curvature tensor diverges. 

The above discussed fixed points and asymptotic behavior have been obtained analytically 
from the equations of motion (20)-(22). Of course it is not possible to find exact expressions 
for solutions which are not at the fixed points. One has to use numerical methods to 
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A 



Y 








y 

Figure 2: Same as figure 1 but for d = 9 



A 



Y 








y 

Figure 3: Same as figure 1 but for d > 9 
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investigate such generic solutions. On figures 1-3 we present some typical results. The 
character of some of the solutions depends on the dimension of the space-time and changes 
at d = 9. 

4 Solutions with finite gravitational constant 

The solutions of the equations of motion discussed in the previous section may be used 
to construct different types of domain wall models. Of course the most interesting are 
those cases for which the lower dimensional effective gravitational constant is finite. The 
effective lower dimensional gravitational constant is obtained by integrating the metric over 
the direction perpendicular to the wall. The are two types of situations when that constant 
is finite. First, the range of the y coordinate may be finite because the space ends at some 
singularities. This type of solutions was discussed for the lowest order string inspired dilaton 
gravity in [6]. Second, the range of y may be infinite but with the warp factor decreasing 
strongly enough to make the integral finite. This type of warped solutions is not present in 
the lowest order theory discussed in [6] but appears in the gravity theory without dilaton at 
the lowest order [5] and in higher orders [4]. 

In the higher order string inspired dilaton theory investigated in the present work both 
mechanisms to get a finite effective gravitational constant can be realized. In fact we can 
define even three classes of such models because there are two different types of singularities 
which can be used to cut off the space. 

Before we present all those classes of solutions a comment on the effective gravitational 
constant is in order. 

In calculating the effective gravitational coupling we assume that the D-dimensional 
metric takes the form 

ds 2 = e 2A ^g { ^dx a dx p + dy 2 . (44) 

and 

= <f>(y) + ip(x a ) . (45) 
We then calculate Z)-dimensional quantities in terms of d dimensional ones: 

D _JAp(<i) 4A a/2 (Jd) n (d) (d) (d)\ ( Af .\ 

K af 3 y s - e K af3y5 - e A [g^g ps - g a5 g M ) ■ (4b J 

and 

Rayfiy = (* 2 + A") . (47) 

Hence we have 3 

Rofi = R { £-e 2A 9 $(dA> 2 + A") 

Ryy = -d(A' 2 +A") 

3 Note that our formulae are different from the ones derived in [12] and used in [10] 
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R = e - 2A R [d) - d ((d + I) A' 2 + 2A") 
R GB = e - 4A R i a ) B -2(d-2)e- 2A R^ ((d - 1)A' 2 + 2A") 

+d(d - 1) (d - 2) ((d + 1) A' 4 + 4A"A' 2 ) (48) 

Then the result of the integration over y of the action (1) gives the string frame d-dimensional 
action 

I d d x^-g( d )e-^ x) (d + C 2 R W + . . .) . (49) 
The d-dimensional cosmological constant is proportional to C\ i.e. to the integral 

J dy exp [dA(y) - <j>(y)] P^A', A", ft, ft'), (50) 

while the inverse of the d-dimensional gravitational constant is proportional to C 2 i.e. to the 
integral 

J dyexp[(d-2)A(y)-<j ) (y)]P 2 (A',A",ft,ft'). (51) 

where P\(A , A" , ft , 0") and P 2 (A' , A" , ft , ft') are some polynomials. If these polynomials 
do not vanish on the solution for other reasons (as happens for the cosmological constant 
on-shell) the behaviour of the exponential allows us to find whether a given solution has 
finite or infinite cosmological and gravitational constants. 

Let us discuss all possible types of solutions which give finite contributions to the effective 
gravitational constant and which can describe the part of the D-dimensional space-time on 
one side of the d-dimensional brane. 

• Type A 

For each d and for each y L there are solutions given by eq. (42) for which \im y ^ yL A' = 
±oo Asymptotically we have 

(d - 2)A' - 0' = . (52) 

y-yo 

Since for d > 2 q < 3/2 for all solutions of eq. (43) therefore the integral (51) is finite 
and these solutions have finite gravitational constant (the corresponding solutions have 
singularities of the type discussed in the lowest order theory in [6]). We can use a part 
of any of those solutions for y G (yz,,yo) to describe the space-time between the brane 
located at y = y and the singularity located at y = yL < yo (the singularity is located 
to the left from the brane). The metric is bounded from above over the whole interval 
{VliVo)- The length of this interval itself is also finite. Thus, the contribution to the 
effective gravitational coupling coming from such solution is finite. 

The above discussed warped metric was obtained by using solutions with A 1 > 
and can describe part of the space-time "to the left" (y < yo) from the brane. The 
analogous metric for (y > y ) can be obtained in the same way from solutions with 
A' < which are related by the symmetries (23), (24). 
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• Type B 



For each d and for each y s there are eight solutions with infinite |A"(y s )| but finite 
values of A'(y s ), four with positive and four with negative A' (they are denoted by 
circles on figures 1-3). For half of them the solution is defined for y > y s and we can 
use part of such a solution for y G (yL = VsiVo) to describe the space-time between 
the singularity and the brane if y s < y . Other solutions are defined for y < y s and 
their part for y G (yo,yn = y s ) can be used when the singularity is to the right from 
the brane. 

All these solutions have singularities at y = y s but of different nature than the solutions 
of type A. Now the metric is non vanishing at y s but the curvature is singular. We 
assume that such a singularity can also be used to cut off the space-time. 



Now we check whether there are solutions without any singularities which nevertheless 
give finite contributions to the effective cosmo logical constant. Any of such solutions 
must be well defined for all values of y > yo (or y < yo). Let us start with solutions 
which approach asymptotically any of the fixed points for y — > oo (y — > — oo). Their 
asymptotic behavior described by eqs. (36)-(40) is enough to determine whether they 
can produce a finite gravitational constant. 

The solution (36) is not appropriate because A' and 0' tend to zero and the warp factor 
approaches a constant for y — > ±oo. The integral over y G (yo, oo) must be infinite. 

The solutions described by eq. (37) also give infinite contributions to the effective 
gravitational constant. We have asymptotically 



After integrating over y we see that this expression diverges when \y\ — > oo for both 
types of solutions: those which approach a fixed point for y — > +oo as well as for those 
approaching a fixed point for y — > — oo. 

Let us now consider the solutions described by eq. (38) for which 



as y — > ±oo. Using eqs. (39) and (41) one can see that the sings of (d — 2) A' — <p' 
and c are the same for 7 = 7+ while they are opposite for 7 = 7_. This means that 
the solutions approaching the fixed points with A' = ±7_/y / a can have finite effective 
gravitational coupling in d dimensions and those with A' = ±^ + /y/a can not. 

After analyzing which solutions approaching asymptotically the fixed points in y = ±00 
lead to a finite effective gravitational constant in d dimensions we have to analyze the 



• Type C 



(d - 2)A' - <f>' - sgn(y) 



a 



(53) 




(54) 
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isolated fixed points. The one with vanishing A' and 4>' (given by eq. (28) for A = 0) 
is not appropriate because it describes just a Minkowski space with a constant metric. 
The situation is more interesting for the two fixed points of the form (28) for negative 
bulk cosmological constant A. The warp factor in the string frame is still constant 
(A' = 0) but the dilaton is not. We get 



J dy exp[(d - 2)A(y) - <f>{y)\ = exp [(d - 2)A - <j> Q ] J dyexp 



a 



(55) 



for (f)' = to with a given by eq. (29). This contribution is finite for the fixed point 
solution with positive 0' = a in the region y G (yo, oo) and for the fixed point solution 
with negative = —a in the region y G (— oo,yo)- Thus the fixed point solution 
A' = 0, 0' = a describes a Randall-Sundrum type warped infinite half-space to the 
right from the brane while that with A' = 0, <f)' = —a the one to the left from the 
brane. 

The situation is similar with the fixed points with non zero A' given by eqs. (30) with 
7 = 7 + . The contribution to the gravitational constant is equal to 



exp [(d - 2)A - O ] J dyexp 




(56) 



Using eqs. (39) and (41) it is easy to check that the sign of the expression in the round 
bracket above is the same as the sign of e^. Thus, the fixed point solutions (30) with 
= +1 (e& = — 1) give finite gravitational coupling after integrating over the region 
y G (yo, +oo) (y G (— oo, yo))- These fixed points are isolated since these conditions on 
€b are just the opposite to the conditions for the existence of asymptotic solutions. 

For 7 = 7_ fixed points with finite gravitational constant are always accompanied 
by asymptotic solutions since (as we have shown) the conditions on e a coincide with 
conditions for the existence of asymptotic solutions. 

We have shown in this section that there are three types of half-spaces which give finite 
contributions to the effective <i-dimensional gravitational coupling. For two of them the 
space-time ends at singularities with divergent curvature. They have different behavior of 
the metric at the singularity, it vanishes for type A while it is non zero for type B solutions. 
Type C solutions have no singularities and extend to infinity. 

All of the above solutions can describe one side of the D-dimensional space-time. Taking 
two solutions we can construct a full space-time on both sides of the brane. There are many 
types of such warped space-times depending on the solutions used to build them. We can 
denote them as: A-0-A, B-0-B or C-0-C where the in X-0-X denotes a brane (which always 
can be put at y — 0) between two solutions of type X (one could consider also non symmetric 
spaces like e.g. A-0-C). The spaces of type A-0-A are similar to those considered in ref. [6] in 
the lowest order theory. Spaces of type B-0-B also end at two singularities but the character 
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of those singularities is different. The curvature becomes singular like in the A-O-A case but 
the metric itself is nonsingular. In case C-O-C the space-time is infinite and warped in such 
a way that the d-dimensional effective gravitational constant is finite (similarly to the space 
of the Randall-Sundrum type). 

Fixed point solutions of the C-O-C type can be written explicitly. Those leading to a 
finite gravitational coupling are given by: 

A' = 0, 0' = sgn(y)4=, (57) 

or 



A = -sgn(y)^ , = -sgn(y) , (58) 



Jol \ \/a 

or 



A = sgn(j/)-= , = sgn(y) , (59) 



a \ \ a 



In all cases the brane is located at y — 0. There are also solutions that asymptotically 
behave as (38) with 7 = 7_ and e a = sgn(y). 

Another interesting type of solutions with finite effective gravitational constant exist for 
all dimensions d different from 9. In figures 1 and 3 one can see solutions which end at two 
singularities of type B. They describe smooth space-times which have finite length along 
the y coordinate without any brane between the singular end points. Such smooth (between 
the end points) brane-less solutions can be described in our notation as B-B. There is also 
a brane-less solutions for d = 9 but it has slightly different character. It can be denoted as 
A-B model because it has a type A singularity at one end a type B singularity at the other 
end. 

One should stress that most of the above solutions appear due to the terms of higher 
order in a'. In the lowest order dilaton gravity only solutions of the type A-O-A are present. 
All other types of solutions (A-B B-B, B-O-B, C-O-C and mixed X-O-Y) appear only after 
the next order corrections are taken into account. 



5 Equations of motion on the brane 

So far we were discussing only the bulk solutions away from the brane. In this section the 
presence of the brane is taken into account. We have to solve the eqs. (15) and (16) including 
Tf£f) and given by (11) and (14). A solution to the left (right) from the brane is denoted 
by A L (y), 4>L{y) (A R (y), 4> R (y)). We assume that A and are continuous across the brane 
while A' and <fi' have jumps described by A' R (y ) — A f L (y ) and <fi' R (yo) — 0l(2/o)- Therefore 
we can integrate these equations over an infinitesimally small interval in y surrounding the 
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brane. Integrating equation a>(3 we get 



e 4>o 



1 - d)A' R + <f>' R -a ( j(d - l)(d - 2)(d - 3)(A' R f - 2(d - l)(d - 2)^0* 



+2(d-l)A^) 2 --(^) 3 



, (60) 



where O — 0(yo)- The dilaton equation gives 



ii 



+ Ct 



2 

-d(d 



l)(d-2)(^) 3 -2d(d 



2/=2/0 



(61) 



Of course (A' R (y ), <f> R (yo)) and (A^(yo), 0l(i/o)) have to independently satisfy the constraint 
equation (21). 

In general it is not difficult to fulfil the above equations for solutions which are not 
exactly at the fixed points. Any of such solutions used to describe part of the space-time 
laying on one side (assume for definiteness y > yo denoted by R side) of the brane has 3 free 
parameters. One of them is the value of A' when approaching the brane A' — > A' R (y ). Of 
course this A' R (yo) cannot be totally arbitrary but can take values only in the range allowed 
for a given solution (see figures 1-3). Specifying the type of a solution and A R (y ) gives 
uniquely a point on a curve at which the brane is located (except for one type of solutions 
when A'(y) is not a monotonic function - then we should specify also the sign of A"(y )). 
The corresponding value of (j)' R (yo) is not a free parameter because it is related to A' R (y ) by 
the algebraic condition (21). 

The remaining 2 free parameters are the constants of integration when we calculate A(y) 
and 4>(y) from A'(y) and <f)'(y), and correspond to the symmetry (23). Those arbitrary shifts 
in A and <fi can be for just one solution reabsorbed in the overall normalization of the metric. 

A pair of solutions used to build a space-time solves the bulk equations of motion for 
1/ ^ !/o and has 6 free parameters: 



A L(yo), A L (y ), (f) L (yo), A' R (y ), A R (y ) , <f> R {y Q ) 



(62) 



The first two requirements, that A(y) and <j)(y) should be continuous at the brane, can 
be used to fix A R (y ) and </)R{yo)- 



A R (y ) = A L (y ) 



(63) 



The remaining two constants of integration, At,(yo) and 0l(|/o), are just the overall normal- 
ization of the metric (but observe that one of them, </> , appears on the left hand sides of the 
brane equations of motion (60), (61)). 

Thus, we are left with two important parameters: A' L (y ) and A' R (y ). These two re- 
maining degrees of freedom are exactly what we need to solve the equations of motion at 
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the brane even for quite general and independent values of the left hand sides of eqs. (60) 
and (61). 4 

The situation is quite different for the solutions which are exactly the fixed point solutions. 
For them A' and 0' are fixed and constant. So a pair of such solution has only 4 free 
parameters 

A L (yo), (/>l(vo), A R (y ), <f> R (y ) . (64) 

Two of them are determined from the continuity of A(y) and (f>(y) across the brane while the 
two remaining can be reabsorbed in the metric normalization. In this case it is clear that 
V((f)(yo)) and V'((f)(yo)) can not be independent because the right hand sides of eqs. (60) and 
(61) depend only on the values of A' and <fi' which are fixed for the fixed point solutions. 

Let us now check what is the relation between V(4>q) and V'(<j)o). Adding eqs. (60) and 
(61) we get 

^ (Vfa) + V'(M) = {A' R [l + 2a ((d - l)(d - 2){A' R f - 2(d - l)A' R 4>' R + («^) 2 )] }|^ 

-{r^l}. (65) 

The expression in the curly bracket above is nothing else as the expression in the curly 
bracket in eq. (25). In section 3 we shown that it vanishes for every fixed point solution. 
Thus the brane equations for fixed point solutions are satisfied only when 

H0o) = -V(<f>o) • (66) 

This means that there are two possibilities to build brane modes using the fixed point 
solutions. One is to use a potential which satisfies eq. (66) for arbitrary value of 0o which 
means that the potential at the brane has the form 

V(<j>) = Ae" (67) 

where A can be interpreted as the brane cosmological constant. The fact that in this case 
the dilaton equation of motion is exactly equal to the potential on the brane (independently 
of the actual value of the dilaton on the brane ) was used in section 2 to argue that with 
such a potential the total effective cosmological constant vanishes. The other possibility to 
solve the equations of motion on the brane is to use more general form of V(<j)) and arrange 
the parameters in such a way that (66) is fulfilled just at = 0o- 

For solutions that asymptotically behave as (38) with 7 = 7. and e a = sgn(y) the r.h.s. 
of (65) does not vanish and can be changed by changing location of the brane. Therefore 
we can arrange it to be equal to the l.h.s. and then the equations of motion are satisfied for 
larger class of brane potentials. 

4 In practice the use of the degrees of freedom (62) can be slightly different, because the l.h.s.s of eqs. 
(60) and (61) do depend on the value of <fi(yo), but number of free parameters is enough to fulfil the brane 
equations of motion for a general potential V{4>). 
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6 Fine-tuning versus (potential) self-tuning of solu- 
tions 



In section 4 we have shown that in the theory defined by the action (1) there are several 
classes of solutions which effectively in ^-dimensions are flat and have finite gravitational 
constant. It is very interesting to check whether those solutions require fine tuning of the 
parameters. 

Let us start with solutions which describe a finite interval in the D-th direction ending 
at some singularities at y — i/l, dr and with a brane located at y — y . The discussion of the 
previous section shows that such solutions do not need any fine tuning at the brane even for 
a general form of the potential V((f>). A pair of solutions used to build a space-time with a 
brane have enough free parameters to fulfil the brane equations of motion (60), (61) for some 
range of values of V((p ) and V'((f>o). This is similar to the so called self-tuning solutions 
discussed for the lowest order gravity with dilaton in ref. [6]. 

The absence of fine tuning at the (central) brane does not unfortunately mean that such 
solutions are fully satisfactory. It was argued in ref. [7] that one should take into account 
also some end point branes necessary to resolve the singularities. Typically the solutions 
have not enough free parameters and the parameters of the end point branes must be fine 
tuned. As a result the whole solution is no longer fine-tuning free. It seems that this kind of 
argumentation can be applied also to the solutions with singularities presented in this work 
for the higher order theory. 

Let us now discuss those solutions which have no singularities (C-O-C in the notation of 
section 4). We have shown in the previous section that the equations of motion at the brane 
(60), (61) can be satisfied for fixed point solutions only when the values of the potential and 
its derivative at the brane are opposite (see eq. (66)) (for special solutions this condition can 
be relaxed). Does this mean that fine tuning is needed for such solutions? The answer to 
this question depends on the definition of fine tuning to be used. Let us look at this problem 
in more detail. 

The problem of fine tuning in the domain wall type solutions in different models has 
been discussed recently in [8, 9]. The fine tuning in these papers was defined as being just 
equivalent to the presence of any relation between values of V and V' at the brane. Some no- 
go theorems have been given in those papers showing that in some theories it is impossible to 
avoid the fine tuning in the brane type solutions with finite effective gravitational coupling 
and without singularities. Our results are in agreement with the spirit 5 of those no-go 
theorems: V((f) ) and V'(0o) are correlated for solutions with finite gravitational constant 
and without singularities. But we do not agree with the interpretation that this means the 
presence of fine tuning. Let us use the following example to explain this. Consider a class 
of theories with the action given by eqs. (l)-(4) with the potential at the brane of the form 
V{4>) = V exp(k(j)). Substituting this to (60) and (61) we can see that the equations of 

5 We write with "spirit" and not with the theorems themselves because none of those papers considered 
the theory discussed in the present paper. 
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motion on the brane can be fulfilled for the solutions of the C-O-C type only for k = — 1. 
According to the definitions from refs. [8, 9] this would mean fine tuning. But k should 
not be treated as a parameter of a given theory but rather as a quantity belonging to its 
definition. Similarly as in the bulk Lagrangian (2)-(3) the same exponential factor exp(— 0) 
multiplies all the terms and we do not treat powers of it as additional free parameters. In the 
case of the bulk Lagrangian this follows form the string theory. It is natural to expect that 
the string theory gives the same factor exp(— 0) also for the brane interactions. Then the 
parameter V can be treated as the brane cosmological constant giving the brane Lagrangian 
in the form 



The condition k — — 1 should not be treated as a fine tuning also for another reason. The 
main feature of usual fine tuning is that a change of one parameter must be compensated 
by changing (fine tuning) some other parameters. It is not the case for the quantity k in the 
above example. For k = — 1 there are solutions of desired properties (for a range of values 
of parameters A and A). But for k ^ — 1 there are no such solutions for any values of other 
parameters. A change of k from —1 can not be compensated by fine tuning of any other 
parameters. 

Taking the above discussion into account we can write the following statement. For 
arbitrary space-dimension D there is a theory defined by the action (1) with the brane 
interactions given by eq. (68) (equivalent to (4) with V{4>) = Aexp(— 0)). In such theory 
there are fixed point solutions which in (D — l)-dimensions are flat, have finite gravitational 
coupling and have no singularities. They are given by eqs. (57) (58) and (59). 

For different brane potentials it is a priori possible to have a solution of the asymptotic 
(38) with 7 = 7_ and e a = sgn(y) but the solution exists only for a window of large negative 
bulk cosmological constant. Those solutions and the fixed point solutions do not need fine 
tuning in the following sense: The free parameters A and A are not correlated. (This is in 
contrast with the Randall-Sundrum model where the bulk and brane cosmological constants 
must be fine tuned against each other.) For a range of values of these parameters there is 
always a solution of this type. 

This way of fulfilling the equations of motion by adjusting some integration constants 
instead of the parameter of the theory are sometimes called self tuning [6] . The better name 
would be rather: potential self tuning. The reason is that the model allows for a self tuning 
(contrary to situation with fine tuning) but the self tuning mechanism itself is still missing. 
We still do not know why this particular type of solutions should be favored by the dynamics 
of the theory. Nevertheless this is a much more interesting situation than that of fine tuning. 

Let us now discuss stability of the (potentially) self tuning solutions with finite gravi- 
tational constant. We have two classes of solutions. One of them consists of isolated fixed 
points (57) and (58). The other consists of solutions that asymptotically approach the fixed 
points (30) with e a = sgn(y) and the fixed point solutions (59) themselves. The isolated fixed 
points are built from the "repulsive" and not "attractive" parts of the fixed point solutions 
in the sense that one cannot perturb the solution and it must stay exactly at the fixed point 




(68) 
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from infinity up to the brane. When such a brane solution of the C-O-C type is slightly 
perturbed it develops a singularity in a finite distance from the brane. This emphasizes even 
more that those solutions are only potentially self tuning ones. 

On the other hand the fixed points (30) with e a — 1 (e a = —1) for y — > oo (y — > — oo) are 
asymptotically approached by a solution (also with finite gravitational constant) that has 
variable A' and 0' - therefore the solution can be perturbed and it would be a much better 
candidate for the self tuning solution if not for the fact that the solution of this type exists 
only for a small window of large negative bulk cosmological constants. 

7 Conclusions 

In the present paper we have investigated the higher order gravity theory coupled to the 
dilaton in arbitrary number of dimensions D. We included in the action the terms with four 
derivatives as motivated by the order a 1 string theory corrections. Equations of motion for 
such theory coupled to a brane have been derived. Their solutions of the domain wall type 
have been analysed in much detail. It occurs that many more types of solutions exist in such 
theory as compared to theories without dilaton or without (some of) higher order terms. In 
particular there are several different possibilities to obtain brane models with finite effective 
(D — l)-dimensional gravitational constant. The space in the directions perpendicular to 
the brane may end with some singularities or may be infinite. There are two types of 
singularities: for one type the metric is non zero at the singularity while for the other it 
vanishes. There is also another new type of solutions which has not been found in previously 
analysed models. These are brane-less solutions with finite effective gravitational constant. 
For those solutions the space ends with two singularities but is smooth between them. 

However the most interesting solutions are those for which the extra direction is infinite 
but the effective gravitational constant is finite due to appropriate warp factors. There 
are to 3 such brane models depending on the value of the bulk cosmological constant 
A. It has been argued [8, 9] that this kind of solutions in similar models always requires 
the parameters of the Lagrangian to be fine tuned. We have shown that in the theory 
considered in the preset work, with all the leading string corrections, it is not necessarily 
the case. Our solutions are fine-tuning free if the bulk potential has a string motivated form 
V(<f)) = Aexp(— 0) where A is the brane cosmological constant. With strong assumptions 
on the bulk cosmological constant (Ao < A < Ai) the solutions are fine-tuning free for 
larger class of brane potentials. The absence of fine tuning means that for arbitrary values 
(within some allowed range) of the bulk and brane cosmological constants, A and A, there 
are solutions with desired properties. They can be obtained by appropriate choice of some 
integration constants. 

The method of adjusting the integration constants instead of adjusting the Lagrangian 
parameters is sometimes called self tuning. But without any mechanism favoring those 
specific integration constants it can be considered only as a potential self tuning. Thus our 
solutions are of the (potential) self tuning type. They exist for a range of the Lagrangian 
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parameters A and A. We discussed also stability of those solutions. It occurs that under small 
perturbations some of them develop singularities in a finite distance from the brane. There 
is however a class of solutions which can be perturbed and does not develop any singularity 
- they however exist only for a small window of large negative bulk cosmo logical constant. 
But anyway to consider them as a solution to the cosmological constant problem we would 
need a mechanism which could explain why out of all possible solutions those which are flat 
and singularity free are favored. Such a mechanism is still missing but those solutions are 
nevertheless very interesting because they at least avoid the fine tuning of the Lagrangian 
parameters. 

One solution similar to those discussed above was reported in ref. [10] when the work 
on the present paper was very advanced. However it is not possible to directly compare our 
results with those presented in [10]. The reason is that the order a' corrections are not the 
same in both cases. We used those derived in the string frame which seems to be the best 
one for considering string corrections. Moreover, the authors of [10] found just one class of 
solutions in a 5-dimensional model while in the present paper more types of solutions are 
presented for arbitrary space-time dimensions. In addition we have discussed stability of 
those solutions and found also many other classes of solutions and their asymptotic behavior. 
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